Anomalous finite size spectrum in the S = 1/2 two dimensional Heisenberg model 
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We study the low energy spectrum of the nearest neighbor Heisenberg model on a square lattice 
as a function of the total spin S. By quantum Monte Carlo simulation we compute this spectrum 
for the s = 1/2, s = 1 and s = 3/2 Ifeisenberg models. We conclude that the nonlinear a model 
prediction for the low energy spectrum is always verified for large enough system size. Ifowever the 
crossover to the correct scaling regime is particularly slow just for the s = 1/2 Ifeisenberg model. 
The possibility to detect this unexpected anomaly with finite temperature experiments on s = 1/2 
isotropic quantum antiferromagnets is also discussed. 
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The square lattice Heisenberg model (HM) has at- 
tracted much attention in recent years because of its 
connection with the antiferromagnetic properties of the 
undoped stoichiometric compounds of High-Tc supercon- 
ductors j^] . The model Hamiltonian in general dimension 
d, reads: 



H 



J ^ Si ■ Sj 



(1) 



<ij> 



where the symbol < i,j > indicates nearest neighbor 
summations, the index i labels the positions Ri of the 
N = L'^ sites on a hypercubic lattice and the quantum 
spin operators satisfy (Si)^ = s(s -I- 1). What makes par- 
ticularly difficult any analytical treatment of this model 
is the fact that the antiferromagnetic order parameter 
TO = jf J2i (with Q ~ {n,TT...)) does not com- 

mute either with the Hamiltonian or with the total spin 
S = Si- Whenever long range order is present in the 
thermodynamic limit, a huge degeneracy of the energy 
spectrum £^(15*) as a function of the total spin S is im- 
plied by the mentioned non commutativity. Hence E{S), 
referenced to the singlet S = ground state energy, is 
predicted to behave as the spectrum of a free quantum 
rotator as long as S" ^ Vn [||-|4| : 



E{S) ^ S{S + 1)/{2IN) 



(2) 



where / is known as the momentum of inertia per site. 

This equation resembles the definition of the spin sus- 
ceptibility X which is however obtained by taking first 
the infinite volume limit of the energy per site e(m) = 
E{S)/N at fixed magnetization m = S/N and then let- 
ting TO 0: e{m) — to^/(2x) + o(to^). An identification 
between / and x is possible only if the excitation spec- 
trum smoothly connects the low energy portion, which 
corresponds to total spin S ^ 0(1), with the regime of 
macroscopic spin excitations: S mN (with to << 1). 
This is a highly non trivial statement which is actually 
verified by the underlying low energy model of the quan- 
tum antiferromagnet (QAF), known as the non linear a 
model (NLctM) §. 



For the HM in one dimension the accepted low energy 
model is instead the Luttinger liquid with K^^ = 1/2 and 
the prediction in this case is that the energy spectrum as 
a function of the total spin behaves as £^(5*) = S"^ /{2xN) 
where x coincides with the spin susceptibility Note 
the strong analogy of the spectrum in one and higher 
dimensions, though in the former case no true long range 
order sets in. 

It is particularly important to verify that, also in pres- 
ence of long range order, the low energy spectrum of the 
microscopic Hamiltonian conforms to the prediction of 
the low energy effective model. This is the main problem 
addressed in this paper where we show the results of high 
accuracy computations of the the energy spectrum E{S) 
in few two dimensional non frustrated spin models. To 
this purpose we use an improved version of the known 
lattice Green's function Monte Carlo (GFMC) technique 
introduced some years ago by Trivedi and Ceperley 0. 
Our method allows to control any form of systematic er- 
ror in a rigorous and simple way For large lattice size, 
the energy difference between spin subspaces becomes ex- 
tremely small. Nevertheless our technique allows a good 
resolution of the spectrum E{S) up to a 16 x 16 lattice 
with a reasonable computational effort. The method we 
use allows to eliminate exactly any source of systematic 
error which affects previous GFMC calculations |^ . For 
large size the results by Runge are also affected by a 
non negligible population control error which is exactly 
eliminated in this GFMC scheme ||]. Our results for 
the ground state energy are instead consistent with those 
shown in Ref. [nOj, obtained with a completely different 
QMC method gj. 

If the quantum top law (|^) is verified, we can calculate 
the unknown quantity 1/2^ on a finite lattice by inserting 
the computed excitation energies at different spin S in 
the equation 



[2xs]-' = N E{S) [5(5+1)]- 



(3) 



Clearly, 1 /xs should approach the physical inverse sus- 
ceptibility for infinite size and for any spin excitation 
S « N provided the quantum top law is verified. 
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FIG. 1. Ax"^ for 2D HM of spin s = 
N=256), s = 1 (data up to N=100) and s : 
N=64). Lines are guides to the eye. 
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It is then instructive to study the deviations from Eq. 
a) through the quantity: 



1 1 



S{S + 1) 



(4) 



which is zero for S = and S — L « N, by definition 
and should vanish for any fixed S in the iV — s- oo limit. 
A finite slope in the numerical data, as clearly shown by 
Fig. |l| for s = 1/2, would imply 1/2/ > 1/2xl, with a 
violation of the quantum top law. A similar and even 
more pronounced, anomaly was also reported in Ref. p^ ] 
to support the existence of a spin liquid phase for the 
spatially anisotropic HM with strong anisotropy. Fig. |l| 
provides clear evidence that this anomaly in the spec- 
trum surprisingly persists in the isotropic 2D HM for 
s = 1/2 HM while is absent, within statistical errors, in 
the higher spin models shown in the same figure. This 
apparently contrasts with the accepted picture stating 
that the low energy properties of the s = 1/2 HM can be 
fully described by the NLcrM. 

How can we reconcile this numerical evidence with the 
established theories of the 2D s = 1/2 quantum antifer- 
romagnets? Does this anomaly have observable conse- 
quence in the low energy behavior of the 2D HM ? In 
order to answer these questions it is important to have a 
better control of the finite size corrections, because the 
data presented in this paper are close to the limits of the 
available computers and considerably larger size cannot 
be handled with the required accuracy. 

In the following, we will make extensive use of spin 
wave theory (^,|l^ (SWT) which is known to be extremely 
accurate for the HM. For instance, the recent quantum 
Monte Carlo (QMC) calculation of the order parame- 
ter gives TO ~ 0.31 p|-p^, quite close to the linear spin 
wave (SW) estimate m ~ 0.3. Moreover SWT, which is 
based on a systematic expansion in 1/s, can be general- 
ized without basic difficulties to finite systems, making 



possible a direct comparison of the numerical results with 
the spin- wave predictions |lj] . 

In order to study the spin excitation spectrum of the 
model as a function of the total spin, we add to the HM 
Hamiltonian H a magnetic field along the z direction: 
Hh = H — hs^Sf . The magnetic field h can be chosen 
to stabilize the desired spin excitation of spin S. In order 
to perform a systematic SW calculation for large s we 
have scaled the magnetic field by s. In this way, when 
s ^ oo the magnetic energy due to the external field 
hs ^ is of the same order ~ as the average energy 
E{S) ^< H >. The classical solution for s ^ oo at fixed 
h is simply obtained by canting the spins of an angle 9 
in the direction of the field (with sin6' = )■ By use 
of the standard Holstein-Primakoff representation of the 
spin operators in terms of bosons and by expanding the 
resulting Hamiltonian in powers of 1/s, we can compute 
the fluctuations over the classical solution and obtain the 
linear SW estimate of the energy of Hh |l^ , |l5|| : 



E{h) = -Js^N{d - h^/8d) - Jds{N - J^k'^k) 



(5) 



^k^h [I6d2(l 
with 7fe 



(cos krr. 



cos k,. 



■■)/d. By expanding 

E{h) about /i = 0, we get E{h) \x^^ {shf which 

immediately gives the well known linear SW correc- 
tion to the classical spin susceptibility xq — I /-id J: 
Xsw/xa = 1 ~ 0.55115/2s [|3|. This calculation can be 
extended to the next leading order and gives the second 
order correction to x- 



x/xo = 1 - 0.55115/2S + 0.0403/(2s)2 



(6) 



The above expansion seems to converge rapidly also for 
s = 1/2. This result differs from the contradicting second 
order SW calculations present in the literature ||l^ . Our 
SW expansion, however, can be in principle different due 
to the scaling of the magnetic field by s, which consider- 
ably simplifies the 1/s expansion of x but has not been 
used before 13 1^. 

The energy spectrum E{S) of the s = 1/2 HM as a 
function of the total spin S can be evaluated within lin- 
ear SWT in two equivalent ways. We can either perform 
the Legendre transform E{S) = E{h) + hsS of Eq. (||) 
at fixed s = 1/2 (referred to as SWH), or we can scale 
the spin S — with fixed /i evaluating the Legendre 
transformed E(fj,s) perturbatively in 1/s before setting 
s = 1/2 (referred to as SWM). Of course these two meth- 
ods are equivalent within a 1/s expansion and the differ- 
ences only arise because we truncate SW expansion at 
a finite order in 1/s. The apparent deviations from the 
quantum top law shown in Fig. can be understood on 
the basis of the aforementioned linear SW analysis. To 
this purpose, in Fig. ^ we compare the SWH prediction 
for l/2xL - 1/2xl/2 with the QMC data. 
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FIG. 2. Size dependence of l/2xs as defined in (^) for 
tlie 2D IfM. Upper curve: linear spin wave. QMC data refers 
to s = 1/2 (circles) s = 1 (squares) and s = 3/2 (triangles). 
Lines are guides to the eye. 

We see that both QMC and SWH results, although 
quantitatively different, show the same dependence on 
lattice size with an apparent plateau at about L ^ 10 
which would suggest a breakdown of the quantum top 
law. While QMC data stop because it is not possible to 
obtain high quality data at larger size, the SW results 
show an abrupt change of behavior with a clear conver- 
gence towards the expected result. The correct 1/L scal- 
ing is however reached only for L > 32. 
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FIG. 3. Ratio of the QMC spectrum with linear SWM 
(lower curve) and SWH (upper curve) for several system sizes. 
The curves are least squares fits of the data. 

The similarity between QMC and SWH, where both 
available, suggests that the non monotonic behavior is a 
genuine feature of the model. In order to provide support 
for this scenario, we considerably reduced the anoma- 
lously large finite size corrections present in Fig. |^ by tak- 
ing the ratio between the energy spectrum evaluated via 
QMC and SWT (using both the SWM and SWH meth- 
ods). As shown in Fig. ^ we get a remarkable collapse of 
all the QMC data for sufficiently large sizes {L ^ 10-^ 16) 



onto a single curve, even at small total spin S. This is 
a direct confirmation that the anomalies detected in the 
spectrum of the s = 1/2 HM are quite similar to those 
present in linear SWT, and eventually disappear only for 
very large size. We stress that this unexpected behavior 
is visible only in the s = 1/2 HM and is strongly reduced 
for higher s (but also in spin- anisotropic systems like the 
XY model). This suggests that quantum fluctuations, 
enhanced by small values of s, deeply affect the proper- 
ties of the model when the spontaneous breaking of a non 
Abelian symmetry occurs. 
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FIG. 4. (a) Ratio between the S=2 and S=l excitation 
energies, (b): Size dependence of l/2xs (^ for S=l (trian- 
gles) and S=2 (squares). The arrow shows the expected l/2x 
(see Fig. The continuous line is a parabolic extrapolation 
(1/2/ = l/2xs for L ^ oo) for 5" = 2 and 6 < L < 16. 

We also notice that at small spin excitations the size 
effects are not so dramatic as shown in Fig. ^(a) where 
the QMC ratio E{S = 2)/E{S = 1) is seen to smoothly 
approach the value 3 predicted by the quantum top law 
(dashed line) . The possibility to control the anomalously 
large finite size effects allows for an accurate numerical es- 
timate of the spin susceptibility even in the s = 1/2 HM. 
In fact, the susceptibility can be directly obtained by ex- 
trapolating to zero magnetization m = S/N the QMC 
results shown in Fig. ||. However, only data representa- 
tive of the thermodynamic limit of the energy per site 
e(m) = E{S)/N must be included in the extrapolation 
procedure. As shown in Fig. ^ the QMC ratios to both 
SWH and SWM predictions can be nicely extrapolated to 
TO = if we exclude points with to < 0.03. On the other 
hand, the thermodynamic limit of the SW results can be 
obtained analytically: e{m) = to^[2xo(1 — 0.55115/2s)]~^ 
and e{m) = TO^pxol^Hl + 0-55115/2s) for the SWH and 
SWM estimate respectively. Both methods consistently 
give the following result for the spin susceptibility of the 
s = 1/2 HM: 



Zx = X/Xo = 0.523 ±0.005 



(7) 



Our value for is also consistent, within error bars, with 



3 



the independent calculation of I shown in Fig. ^(b) and is 
not far from the second order S W value ^ . This results 
slightly differs from series expansion Ref. but much 
more from previous finite size scaling analysis based on 
QMC data on similar lattice sizes ppO| . 
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FIG. 5. Same as Fig. ^(b) for the spin excitation S — L. 
Curves are parabolic extrapolations. The arrows are the ex- 
pected infinite size limits. 

The latter discrepancy can be easily understood on the 
basis of the strong finite size anomaly in the spectrum of 
the s = 1/2 HM we have discussed. In fact, as shown in 
Fig. ^(a) if we extrapolate to the thermodynamic limit 
the known SW spectrum in lattices with L < 16 by use 
of the correct asymptotic scaling xl = Xoo + a/L + h/L? 
[Q,^, we get a poor estimate of the exact L ^ oo limit 
of Xsw- In fact, SWT suggests that the correct scahng 
is valid only for extremely large sizes (L > 64) while it is 
not accurate for the lattice sizes where QMC simulations 
are available (see Fig. ^(b)). 

A further evidence that the data shown in Fig. ^ are 
very weakly size dependent comes from the consistency 
with the NLcrM next leading m ^ prediction Q for 
e(m) = (l/2x)m2 - {l/UTrc^x^)'!^^ ■ This allows to de- 
termine from the slopes for m ^ in Fig. || the SW 
velocity c/cq — 1.18(4) (upper curve), c/cq — 1.1(10) 
^>wer curve) in good agreement with the expectations 

In conclusion we have given convincing evidence that 
the low energy spectrum of the 2D quantum antiferro- 
magnet is consistent with the prediction of the (2+l)D- 
NLcrM. Unexpected anomalies in the finite size spectrum 
are however detected in the s = 1/2 HM. The physical 
origin of these size corrections can be understood on the 
basis of linear SW approximation: One of the two gap- 
less excitations present in the SU(2) model, acquires a 
mass when an external magnetic field is applied. The 
competition between the massive and the massless mode 
becomes quite strong in the m limit leading to the 
observed anomalies. Making use of the results of SWT, 
we have been able to obtain a very accurate estimate for 



the spin susceptibility. 

This extensive discussion of finite size effects in the low 
energy spectrum of the 2D-HM has interesting connec- 
tions with the finite temperature behavior of the system. 
In fact, on the basis of the mapping to the NLctM, the 
properties of the HM at finite temperature T coincide 
with those at T = of a HM in an infinitely long strip 
of width L = 1/T and periodic boundary conditions. By 
applying the finite size SW analysis to the stripe geome- 
try, we then evaluate the uniform magnetization m in an 
external magnetic field h at temperature T: 



m{h, T) ~ xh + 



hT 



ln(l 



-h/T^ 



(8) 



This expression, asymptotically valid for small h and T 
but arbitrary ratio hlT ^ obeys the known scaling form 
predicted by Fisher The anomalous finite size cor- 
rections in the energy spectrum we have uncovered also 
suggest the possible presence of a small crossover tem- 
perature Tx corresponding to the maximum in Fig. ^. In 
the HM, ksTx can be estimated as few percents of the 
natural energy scale J. The observable consequences of 
this possible anomalous behavior of the low temperature 
susceptibility of the 2D HM should be reconsidered on 
the basis of this study. 
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